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Abstract 

By numerically solving the quantum kinetic equations we compute the range 
of parameters where the v„ — > u s oscillation solution to the atmospheric neu- 
trino anomaly is consistent with a stringent big bang nucleosynthesis (BBN) 
bound of Nffi N ~ 3.6. We show that this requires tau neutrino masses in the 

range m„ T ~ 4 eV (for |5m| tm | = 10~ 2 ' 5 eV 2 ). We discuss the implications 
of this scenario for hot+cold dark matter, BBN, and the anisotropy of the 
cosmic microwave background. 
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I Introduction 



Neutrino physics continues to provide the most promising window on physics beyond 
the standard model. There are numerous indications for neutrino oscillations including 
the atmospheric neutrino anomaly, solar neutrino problem, LSND experiment, structure 
formation in the early Universe, indications of NfM N < 3 etc. For a recent review, see e.g. 
Ref. [0. 

The atmospheric neutrino anomaly P] has been confirmed by the superKamiokande 
experiment || . The observed up-down asymmetries of the detected muons indicate Q that 
the simplest solution to this anomaly is either — * v T |5|] or z/ M — > u s [|5|-|T(]|] oscillations 
(although significant additional mixing with v e cannot currently be excluded if Sm 2 ~ 2 x 
10~ 3 eV 2 [|11|]). In each case the oscillations are maximal or nearly maximal (0.80 ~ sin 2 26 ~ 
1) and 

4xl0- 4 (10- 3 )^|5m 2 tmos |/eV 2 ^10- 2 , (1) 

for — > v T (z/ M — > v s ) oscillations Qj. Although experimentally similar, the — > v T and 
v s oscillation solutions to the atmospheric neutrino anomaly can be experimentally 
distinguished (l^Q- These two solutions will also have quite different implications for early 
Universe cosmology. Naively the — > v s oscillation solution to the atmospheric neutrino 
anomaly would appear to lead to an effective number of four neutrinos [[3j|] in the early 
Universe, which would make it difficult to reconcile big bang nucleosynthesis (BBN) with 
the estimated primordial light element abundances (for a recent discussion, see Ref. [|R|]). It 
should be stressed that the ordinary-sterile neutrino oscillations can only populate the ster- 
ile state provided that the lepton number of the Universe is very small. Since the origin of 
the baryon and lepton asymmetries of the Universe are unknown it is possible that a lepton 
asymmetry was created at some early time T 3> 100 MeV . Furthermore, if the lepton to 
photon ratio is larger than about 10~ 5 then the BBN bounds will be evaded [0. Clearly 
this is one possibility. However another possibility is that the neutrino oscillations them- 
selves generate the lepton number. In fact, a careful study of the ordinary-sterile neutrino 
oscillations in the early Universe reveal that the oscillations themselves typically generate 
a much larger asymmetry for a large range of parameters [|18| -|2~THPI- As already discussed in 



Ref. [2~0|, there is a concrete mechanism whereby the large neutrino oscillation generated 
asymmetry can prevent the (near) maximal — ► v s oscillations from populating the sterile 
neutrino in the early Universe. This scenario requires the tau neutrino to be in the eV mass 



x There is an astrophysical/cosmological arguement which favours Un — ► v s over v„ — ► v T . Ref. [13] 
points out that the decaying neutrino theory for the ionisation of hydrogen in the interstellar 
medium, in conjunction with the assumption that the cosmological constant is zero favours — ► v a 
over — ► v T . 

2 It turns out that this asymmetry is typically generated at temperatures less than 100 MeV 



so that the quantitative calculations of Ref. [17| are not generally applicable to the case where 
neutrino oscillations generate the asymmetry. 
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range and to mix slightly with the sterile neutrino. In this case the v T — > v s oscillations 
generate a large tau neutrino asymmetry in the early Universe which suppresses — > v s 
oscillations (the tau neutrino asymmetry can suppress — > v s oscillations because the 
matter term for — *> v s oscillations also contains a part which is proportional to the tau 
neutrino asymmetry of the background plasma). The calculation is somewhat non-trivial 
because the — > v s oscillations like to produce a large muon neutrino asymmetry which 
can compensate for the effect of the large tau neutrino asymmetry in the matter term for 

— > v s oscillations. The computation of Ref. [fJ(J utilised an approximate solution to the 
quantum kinetic equations (which was called the 'static approximation' in Ref. PU| ) and 
is valid provided that the system was smooth enough. This is a very useful approximation 
because it saves considerable CPU time as well as giving more insight than the more com- 
plicated quantum kinetic equations. However as pointed out in Ref. |20[ this approximation 
is not always valid for the entire range of sin 2 2# , Sm 2 parameter space of interest. In 
particular, if sin 2 29 is large enough then this approximation is generally not valid because 
the lepton asymmetry is created so rapidly. Thus, the main purpose of this paper is to do 
a more exact computation (by using the quantum kinetic equations rather than the simpler 
static approximation) of the region of parameter space where the — ► v s solution to the 
atmospheric neutrino anomaly is consistent with BBN (assuming Nf^ N ~ 3.6). Also we 
will also include the entire range Eq.([l|) of 5m 2 atmos . If the sterile neutrinos are not pop- 
ulated by — > v s oscillations then during the low temperature evolution of the neutrino 
asymmetries some electron neutrino asymmetry will be transferred from the tau neutrino 
asymmetry (due to v T — > v e oscillations), which significantly affects BBN. This observation 
was first made in Ref. |2TJ and we include a discussion of this effect here for completeness. 
We will also discuss the implications for the hot+cold dark matter model and the anisotropy 
of the cosmic microwave background (CMB). 

For the purposes of this paper, we assume that there is only one light sterile neutrino. 
Thus, we are considering a four neutrino model. Probably the simplest four neutrino model 
which can explain the atmospheric neutrino anomaly through — > v s oscillations, and also 
explain the solar neutrino problem is the model of Ref. [[| . In this model the solar neutrino 
problem is explained by the oscillations of v e —>■ u^, v s with parameters consistent with the 
small angle MSW solution of the solar neutrino problem. Our results (as well as the previous 
work pOl^H) is applicable to this model. However this model does not seem to be compatible 
with the LSND anomaly [^2|. Actually, as we will discuss in more detail in section VIII, our 



results indicate that Nfj^ < 4 does not seem to be consistent with any four neutrino model 
which explains all three experimental anomalies. Hence, if experimental data indicate that 

— > v s oscillations are required to explain the atmospheric neutrino data, and if the solar 
and LSND anomalies have been correctly interpreted in terms of neutrino oscillations, then 
N^ff < 4 actually suggests the need for more than four neutrinos. 

Admittedly, the theory of neutrino oscillations in the early Universe is quite a complicated 
subject. The readers who are primarily interested in our results can skip directly to Figures 
2,3,4. These figures summarise the main results of this paper. 

The outline of this paper is as follows: In section II we briefly review the phenomenon of 
neutrino oscillation generated neutrino asymmetry in the early Universe. In section III we 
explicitly write down the quantum kinetic equations which we will need latter on, and we 
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numerically solve them for some illustrative examples. In section IV we obtain the region of 
parameter space where the — > v s oscillation solution of the atmospheric neutrino anomaly 
is consistent with a BBN bound of N^ N ~ 3.6. In section V, VI, VII we discuss the detailed 
implications of this four neutrino scheme for BBN, hot+cold dark matter model, and the 
anisotropy of the CMB. In section VIII we conclude. 

II Oscillation generated neutrino asymmetry in the early Universe 

Our notation/convention for ordinary-sterile neutrino two state mixing is as follows, the 
weak eigenstates (u a , u s ), with a = e,/iorr, are linear combinations of two mass eigenstates 

v a = cos 9 u a + sin 9 u^ u s = — sin 9 Q v a + cos 9 u b . (2) 

Note that we define 9q such that cos 29q > and we take the convention that 5m 2 = m 2 —m 2 . 
The neutrino asymmetries are defined by, 

L s r^n^ (3) 

with rii being the number density of species i. Finally, note that when we refer to a neutrino, 
sometimes we will mean neutrinos and anti-neutrinos. We hope the correct meaning will be 
clear from context. 

It was shown in Ref. that for v a — > v s (a — e,fi, or r) oscillations (with \Sm 2 \ ~ 



10 eV ) the evolution of lepton number has the form (for small L Va ) 



dL -~c(L„„ + ?), (4) 



dt \ a 2 

where r] is set by the relic nucleon number densities (and is expected to be small, rj/2 ~ 10~ 10 ) 
and C is a function of time, t (or equivalently temperature, T). At high temperature C is 
negative, so that (L Ua + rj/2) ~ is an approximate fixed point. However if 5m 2 < 0, then 
C changes sign at a particular temperature T = T c . At this temperature rapid exponential 
growth of neutrino asymmetry occurs (unless sin 2 29 is very tiny, see Eq.(H) below). The 



temperature where C changes sign was calculated to be [18 



i 



/-5m 2 cos2#oV T , 
T c ^16(^ ^y^J MeV. (5) 

The generation of neutrino asymmetry occurs because the v a — ► v s oscillation probability is 
different to the v a — > v s oscillation probability due to the matter effects in a CP asymmetric 
background. As the asymmetry is created, the background becomes more CP asymmetric 
because the neutrino asymmetries contribute to the CP asymmetry of the background. 

The phenomenon of neutrino oscillation generated neutrino asymmetry was studied in 
more detail along with some applications in Refs. p9j- f2T| . |23| -[25[] . Ref. [B5 calculated the 



region of parameter space where the neutrino asymmetries were created under the approx- 
imation that all the neutrinos have a common momentum (p ~ 3.15T). In Ref. [2D|, the 



4 



neutrino momentum spread was taken into account by developing an approximate solution 
to the quantum kinetic equations. Further studies of this approximation, along with an 
alternative derivation appears in Ref. p3 |. Ref. [J2(J also contains a detailed study of the 
temperature region T ~ T c where the initial exponential growth of the neutrino asymmetry 
occurs. In Ref. []ST| the low temperature T ~ T c evolution of the neutrino asymmetry was 
studied. After the initial exponential growth of the neutrino asymmetry occurs, the colli- 
sions and eventually MSW transitions combine to keep the asymmetry growing. This low 
temperature evolution of the asymmetry is approximately independent of sin 2 29 (assuming 
that sin 2 29 <C 1). The 'final' value of the asymmetry arises at the temperature ]2T] 

Tt ~ 0.5(\5m 2 \/eV 2 ) 1/4 MeV. (6) 



The magnitude of the final value was calculated to be |2l 

L f v /h ~ 0.29 for \5m 2 \/eV 2 ~ 1000, 



Ll/h 
Ll/h 



0.23 for 3 ~ \5m 2 \/eV 



2 < 



1000, 



0.35 for 10~ 4 ~ \5m 2 \/eV 2 ~ 3, 



(7) 



where h = (T 3 /T 3 ). Strictly these results are valid for v a — * v s oscillations in isolation (i.e. 
in the idealized case where there are only the two flavour oscillations v a — > v s occurring). For 
the realistic case of three ordinary and one sterile neutrino considered in this paper, these 



results hold approximately for the v T — > v s oscillations (assuming that m 2 T 3> m? h 
This large neutrino asymmetry occurs for a wide range of sin 2 29 , 5m 2 P0,[l9~ 



,m: 



5 x 10 



-10 



eV 2 
1 5m? 



~ sin 2 ' 



2# ~ 4(2) x 10" 



eV 2 
1 5m? 



\5m 



> 



10" 4 e^ 
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(8) 



with 5m? < for v^^ — > v s [y e — > z/ s ) oscillations. Note that the upper bound on sin 2 29q 
in Eq.(^|) assumes that the energy density of sterile neutrinos, which arise in the period 
before the exponential growth of neutrino asymmetry occurs (i.e. T ~ T c ) is less than 0.6 
of a standard neutrino species. That is 6N*/j N ~ 0.6 is assumed in Eq.(|8|). Because of the 
present observational uncertainties in the primordial light element abundances this bound 
is not rigorous. It is nevertheless interesting to suppose that the bound is rigorous and 
to explore the consequences. In any case it is clear that the generation of large neutrino 
asymmetries is quite a general phenomenon which occurs for a large region of parameter 
space if light sterile neutrinos existfj]. 



3 The alternative case where and v T are approximately maximal mixtures of two mass eigen- 
states is discussed in Ref. [24]. 

4 In the region of parameter space where |<5m 2 | <ti 10~ 4 eV 2 , the evolution of lepton number 
is dominated by oscillations between collisions and the lepton number tends to be oscillatory 
[26,27,il|. Note however that, if say v T — > v s oscillations have \5m?\ ~ eV 2 , then the final 



value of L Vt is so large that oscillations of say v e — ► v s with \5m | -C 10 eV are heavily 



suppressed by the large matter effects (caused by the large L„ T ) for T ~ 0.5 MeV, and thus such 
oscillations cannot have any effect for BBN. 
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Oscillation generated large neutrino asymmetries can have several important implications 
for cosmology, including, 



(1) The BBN bounds on ordinary - sterile neutrino oscillations can be evaded p0|p3 



(2) If electron asymmetry is generated, then big bang nucleosynthesis can be directly affected 
though the modification of nuclear reaction rates (such as v e + N <-» P+e~ , v e + P <-» iV+e + ) 
2T]& 



(3) The modification of the neutrino number densities will have implications for structure 
formation and will affect the anisotropy of the cosmic microwave background. 

In this paper we will illustrate all three of these implications within the framework of the 
four neutrino model which solves the atmospheric neutrino anomaly by — > v s oscillations. 

Ill The quantum kinetic equations for ordinary - sterile neutrino oscillations in 
the early Universe 



The density matrix [p8-30| for an ordinary neutrino, v a (a = e, /x, r), of momentum with 
magnitude p oscillating with a sterile neutrino in the early Universe can be parameterised 
as follows, 

PvM = \Po(p)[I + P(p).*], PM = \hip)[I + P(p).a], (9) 

where I is the 2x2 identity matrix and P(p) = P x (p)^t+ P y (p)y + P z (p)z, a = a x x+a y y+a z z 
(the Oi are the Pauli matrices). It will be understood that the density matrices and the 
quantities Pj also depend on time t or, equivalently, temperature T (the time temperature 
relation for m e ~ T ~ m M is dt/dT ~ -M P /5.5T 3 , where M P ~ 1.22 x 10 22 MeV is the 
Planck mass). 

We will normalise the density matrix such that the momentum distributions of u a , v s 
are given by 

N Va = \Po(p)[I + Pz(p)]N°(p), N Vs = \Pv(p)[I - P z (p)]N°(p), (10) 

where 

N°(p) = ^— r-r- (11) 

^ 27T 2 1 + 6Xp (|) 1 J 

Note that there are analogous equations for the anti- neutrinos (with P(p) — > P{p) and 
Po —* Po)- The evolution of Pq(p) and P(p) [or Pq(p), P(p)] are governed by the equations 



= V(p) x P(p) + [1 - P,(p)][^nP (p)]z 
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[D{p) + ^lnP Q {p)][P x {p)x + P y (p)y], (12) 



K(p)-±P (p)(l + P z (p)) 



(13) 



where D(p) = T(p)/2 and T(p) is the total collision rate of the weak eigenstate neutrino of 
momentum p with the background plasmaQ and K(p) = N eq (p) / N° (p) , with N eq (p) being 
the expected number of neutrinos in thermal equilibrium, i.e. 

N eq ( P ) - ~ 4^r-v (14) 



2tt 2 i + exp [2±^ ■ 



For anti- neutrinos, \i Va — > /ip a in the above equation. The chemical potentials /jL Va , fip a , 
depend on the lepton asymmetry. In general, for a distribution in thermal equilibrium 

T^rlnr 1 poo T^rlnr 

(15) 



4£(3) Jo l 4- e x+va 4£(3) Jo l 4- e x+iZ* ' 
where /i a = ^ Ua /T and /z s = /i Pc /T. Expanding out the above equation 



L ~ 



24C(3) 



(16) 



which is an exact equation for /^ a = — /^, otherwise it holds to a good approximation 
provided that 1. For T ~ T£ c (where T| ec ss 2.5 MeV and T^ T C ss 3.5 MeV are the 

chemical decoupling temperatures), ~ — Aip a because processes such as v a + v a e + + e~ 
are rapid enough to make ^ a + ~ /^ e + + ^ e -— 0. However, for lMeV ~ T ~ TJ^ C , 
weak interactions are rapid enough to approximately thermalise the neutrino momentum 
distributions, but not rapid enough to keep the neutrinos in chemical equilibrium (]. In this 
case, the value of H a is approximately frozen at T ~ T£ ec (taking for definiteness L Va > 0), 
while the anti-neutrino chemical potential H & continues increasing until T ~ 1 MeV. 

The form of the evolution equation, Eq.(|T2]) is very easy to understand. The VxP term is 
simply the quantum mechanical coherent evolution of the states (this is derived in many text 
books, see e.g. Ref. P5| ). The damping term (D[P x St + P y y]) is just the destruction of the 
coherence of the ensemble due to collisions. The only suprising thing is that D(p) is half of the 
collision frequency i.e. D(p) = T(p)/2 [rather than T(p)] [02]| . The rate of change of Pq is due 
only to the repopulation of ordinary neutrinos due to weak interactions (since, from Eq. (|T0|), 
Po{p) = (N Ua + N Us )/N°, which is obviously unchanged by the v a — > u s oscillations). Note 



5 From Ref. [HJH it is given by T(p) = yG 2 F T 5 (p/3.l5T) where y ~ 4.0 for v = v e and y ~ 2.9 
for v = Vu, v T (for m e ~ T ~ mj. 

6 The chemical and thermal decoupling temperatures are so different because the inelastic collision 
rates are much less than the elastic collision rates. See e.g. Ref. for a list of the collision rates. 
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that Eg. (p~3|) is approximate and holds provided that the distributions of the neutrinos and 
background particles are in near thermal equilibrium (see Ref. [25| for a detailed derivation of 
this formula). Finally the dlnP /dt terms arise because the repopulation does not populate 
sterile states or mixtures of states. [For example, because the repopulation does not affect 
the number of sterile states, it follows that ^jf\ r e P o P uiation = (1 — Pz) 91 q^° , since |-Po(l — P z ) 
is unchanged by the re-population]. 

The quantity V(p) is given by [P9LB0] 



V(p)=/3(p)x + A(p)z, 



where f3{p) and \{p) are 



sin 29 , A(p) 



8m 2 
~2p~ 



[cos29 -b{p)±a(p)}, 



(17) 



(18) 



in which the +(— ) sign corresponds to neutrino (anti- neutrino) oscillations. The dimen- 
sionless variables a(p) and b(p) contain the matter effects (more precisely they are the 
matter potential divided by 5m 2 /2p). For v a — > v s oscillations a(p),b(p) are given by 



a(p) 



-4((3)V2G F T 3 L^p ^ -4C(3) v / 2G F TM a p 2 



7T 



l 8m 2 



n 2 5m 2 M? 



(19) 



w 



where £(3) ~ 1.202 is the Riemann zeta function of 3, Gf is the Fermi constant, Mw is the 
W— boson mass, A e ~ 17 and A^ T ~ 4.9 (for m e ~ T ~ m M ). The quantity is given by 

=L Va +L Ve +L^+L UT +m (20) 

where 77 is a small term due to the asymmetry of the electrons and nucleons and is expected 
to be very small, rj ~ 5 x 10~ 10 . Recall that the neutrino asymmetry is defined in Eq.(^). 
For example, the number density of v a is 



poo roc 1 

/ iV, Q dp= / -P (l + P z )iV>. 
jo jo 2, 



(21) 



Note that in the interests of notational simplicity, in the above equation and in the following 
discussion, the functional dependence of Pi, N°, (3, A, D on the neutrino momentum will not 
always be made explicit. The rate of change of lepton number is given by 



dL Va 


d 


n v — n p 


d 


n Vs - n Us 


dt 


~ dt 


n 7 


~ ~dt 


n 7 



(22) 



Thus, using Eq . (|10|) 



dL h 



dt 



2n„ 



dPo 
dt 



1-P* 



OP 1 _dPo 

dt dt 



i-p, 



p dP : 



dt 



N°dp. 



(23) 



It is possible to numerically integrate the system of coupled differential equations, 
Eq. (|12 , |13 ,23), although it is (CPU) time consuming. Unfortunately we do not have a par- 
allel supercomputer handy. We therefore employ the useful time saving approximation of 
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integrating the oscillation equations, Eq.(|T^), in the region around the MSW resonance (we 
will define this region precisely in a moment). Actually away from the resonance the oscil- 
lations are typically suppressed by the matter effects or by sin 2 29 (or both). Thus, this 
should be a good approximation, which we can check by taking larger slices of momentum 
space around the resonance. 

In the adiabatic approximation, the oscillation probability for v a —>■ u s is given by the 
formula, 

P = sin 2 2# m (sin 2 -^->. (24) 

Of course this formula is only valid provided that sin 28 m is approximately constant on the 
interaction time scale (1/D). In Eq.(^), the brackets (...) denotes the average over the 
interaction times r, and sm 2 29 m ,L m are the matter mixing angle and oscillation length 
respectively. In terms of (3, A, D, 

• 2 9/1 _ P t2 _ 1 

8111 Wm ~ 02 + A 2> L m ~ 02 + A 2> 

,2 ' \ _ n / „-tD -2 ' 



(sin 2 ) = D / e~ rD sin 2 dr. (25) 

2 1~ J TYl 2 T J wyy 

It is straightforward to show that [^(J 

1 B 2 

V = t\ (26) 

2 3 2 + D 2 + \ 2 ' 1 ; 

Thus the oscillation probability has the MSW resonance when X(p) = (note that for 
notational convenience we call this the 'MSW resonance even if D 2 3> f3 2 ). Solving the 
equation, \{p = p res ) = we find, 



X 



2 



Pres r> ~^ 



( X 2 \ 2 cos 2^ 



+ (27) 



2X 1 \| \2Xj X x 

where Xx = b(p)/p 2 and X 2 = a(p)/p (note that Xi are independent of p). The resonance 
width in momentum space, A, can be obtained from X(p = p res ± |A) 2 w /3 2 + .D 2 , 

A * 4p ™ (28) 
|5m 2 |(2p res X 1 -X 2 )' 1 ; 

where /3, D are evaluated at p = p res . For anti-neutrinos, X 2 — > — X 2 in Eqs. (^7|,|28]) . 



For the numerical work the continuous variable p/T is replaced by a finite set of momenta 
x n = p n /T (with n = 1, 2, N) on a logarithmically spaced mesh. The variables P(p) and 
Pq(p) are replaced by the set of N variables P(x n ) and P (x n ). The evolution of each of 
these variables is governed by Eqs.(|l^,|l3|), where for each value of n, the variables V(p) 
and D(p) are replaced by \(x n ) and D(x n ). The anti-neutrinos are similarly treated. Thus 
the oscillations of the neutrinos and anti-neutrinos can be described by 8N simultaneous 
differential equations. For each time step the lepton number is computed from Eq.(|23|) and 
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the chemical potentials are obtained from Eg. (|Tq) . In order to integrate first order differential 
equations initial conditions must be specified. At very high temperatures the oscillations 
are heavily suppressed by the collisions so we take P x>y {p) = 0,P Xty (p) = initially [to see 
this observe that at high temperature, D 3> (3 so that V — > from Eq. (p6|) 1. Also, we 
assume that there is initially a negligible number of sterile neutrinos which means that we 
take P z (p) = P z {p) = 1 and Po(p) = Po{p) = 1 initially. The value of our initial temperature 
must be significantly higher than the temperature where the exponential growth of lepton 
number occurs [we typically used T initia i ~ 4T C , where T c is given in Eq.(|5])]. In our numerical 
work we integrate Eq.(|T2|) around the region 

A A 

Pres - f— < P < Pres + f~, (29) 



with f = 7. Note that repopulation equation Eq.(|l"3"D needs to be integrated over the entire 
region in momentum space, which we typically approximated to be 0.001 < p/T < 20. We 
also take into account the damping in the momentum region away from resonance. That is 
for the momentum region, 

A A 

P < p res - f— and p > Pres + f — , (30) 

we neglect the oscillations and Eqs. ( jT2|Ji3D are truncated to, 

^^ = -(m+^inP (p)y x M, 
(jR{p) -(i-p,(p))JWp), 



dt y ^" dt 

(31) 



Dt 



K(p)-±P (p)(l + P z (p)) 



and similarly for P Xj y >z (p), Pq{p). We checked the stability of our results by integrating 
Eq. (|12"D over larger slices of momentum space [i.e. taking / > 7]. 

We have numerically solved these equations for three illustrative examples. In Figure 1 
we have plotted \L„ T \/h versus T (recall that h = T^/T^) for v T — > v s oscillations. We have 
chosen Sm 2 /eV 2 = —0.5, —50, —5000 and sin 2 29o = 10~ 8 with initial condition L Va = 
at TinUiai = 300 MeV. In these examples the rapid exponential growth of lepton number 
occurs when T c ~ 16, 37, 80 MeV. In the region T > T c , Z>) -> is an approximate fixed 
point and this explains why L Ua — > —rj/2 ~ —2.5 x 10 -10 . In fact this behaviour occurs 
provided that the initial value of L Va is less [] than about 10~ 5 . 



7 For initial values of L Ua ~ 10~ 5 the oscillations are not able to drive — > 0. In this case the 
oscillations do not significantly destroy the initial asymmetry in the region T > T c . In the lower 
temperature region T < T c the oscillations increase the size of the asymmetry and the final value 
of the asymmetry is also given approximately by Eq.(|7]) (assuming of course that the magnitude 
of the initial asymmetry was less than L[ to begin with). 
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The behaviour illustrated in Figure 1 is well understood 0,^1 an d quite general provided 
that 5m 2 < 0, \5m 2 \ ~ 10~ 4 eV 2 and sin 2 26$ is within the wide range, Eq.(^). In the region 
before the exponential growth of lepton number occurs (T ~ T c ), I>) « (which means 
that X 2 ~ 0) so that the MSW oscillation resonance for neutrinos p v res has approximately 
the same value as the MSW oscillation resonance for anti- neutrinos (f£ es )- As the neutrino 
asymmetry is created (taking L Ua > for defmiteness) , p^ es /T ^> 1 (and typically ~ 20) 
and f r jT < 1 (and typically ~ 0.2)0. In the region after the exponential growth (T ~ T c ) 
the collisions keep p v res jT < 1. For lower temperatures T ~ T c /2, the b(p) term [in Eq. (|T^)1 
can be approximately neglected (relative to the a(p) term) and 

v v 1 1 

£Z££ ~ n _ (32) 

T ~ TX 9 T 4 LW ' 1 ' 



As T increases, p v res jT increases and MSW transitions convert v a — > v s at the momentum 
p = p v res . The effect of this is to keep L Va growing until all of the v have passed through 
the MSW resonance^]. If there were no re-population effects (which is approximately true 
for the case of low \5m 2 \ ~ 3 eV 2 , where the asymmetry doesn't become large until low 
temperatures, T ~ 1 MeV) then the expected asymmetry is L Ua ^ n„ a /n 7 ~ 0.375 (since 
all of the v a have been converted into v s ). However in the case where L Va becomes very large 
in the region above about 1 MeV, repopulation effects must be taken into account, and the 
final value of the asymmetry is typically reduced somewhat [see Eq.(^)]. Thus, this explains 
why the final asymmetry is so large and approximately independent of sin 2 26$ so long as 
it is small. Note that in the case of relatively large sin 2 29q, significant number of u s , v s are 
populated in the region T ~ T c . In this case, the MSW transitions at low temperature are 
less effective at creating L Ua . Indeed, in the limit where the u s , v s are fully populated, MSW 
transitions would simply interchange equal numbers of v a with v s and thus there would be 
no significant generation of neutrino asymmetry in this case. 

Note that for the examples in Figure 1, the sign of L Vr changed at the temperature 
T ~ T c . This behaviour is expected (see Ref. ]2(| for a discussion of this point). Actually our 



0) 



numerical integration of the quantum kinetic equations reveals that for very large sin 2$, 
oscillations of L Ut occur. In this case, it would presumably be impossible to predict the sign 
of L Vt (see also the discussion in Ref. |l|§). Finally note that this result was not found in 



Ref. |^0[ because the static approximation developed there is not valid for sin 2 26$ sufficiently 
large (due to the extremely rapid exponential growth of lepton number). 



8 See Ref. for a figure illustrating the evolution of p^es f° r an example. 

9 Note that the MSW transitions keep p^ es /T approximately constant from T ~ T c /2 until L Ut 
is quite large (~ 10~ 2 ). From Eq.(^2|), p v res jT m constant, implies that L Ut is proportional to 
1/T 4 . This explains why -g^j-^p- — —4, in the region after the exponential growth of lepton number 
occurs for the examples in Figure 1. 
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IV Consistency of — > u s solution of the atmospheric neutrino anomaly with 
BBN 



We now turn to the main issue of this paper, that is, to accurately determine the region of 
parameter space where the — > v s oscillation solution to the atmospheric neutrino anomaly 
is consistent with a stringent BBN bound of N^M N ~ 3.6. 

Let us begin by introducing some notation. We will denote the parameters, sin 2 2#o, Sm 2 
and the matter terms a(p),b(p) [Eq.(|iTS|)] for v a — > v s oscillations by sin 2 29q s , 5m 2 as and 
a as (p) , b as (p) respectively (a = e, /x,r). We are assuming that there is only one sterile 
neutrino which mixes maximally or near maximally with the muon neutrino so that the 
atmospheric neutrino anomaly is solved, i.e. 

sin 2 2C ~ 1, 1(T 3 - \Sml s \/eV 2 - 1(T 2 . (33) 

We assume that 



m Ve , m Vti , m Ua < m„ T , (34) 

and that the v T oscillates with the v s with small mixing, i.e. sin 2 29™ <C 1. In fact it turns 
out that the above assumptions are actually necessary if NfM = 4 is to be avoided. 

Of course there will be many specific particle physics models consistent with the assump- 
tions, Eqs. (|33|j34] ). For example, Ref. [[J discusses one such model where in addition to ap- 
proximately maximal — > v s oscillations there are v e — > z/ M , v s oscillations with parameters 
consistent with the small angle MSW solution of the solar neutrino problem (sin 2 26*o ~ 10 -2 , 
\5m 2 \ ~ 10~ 5 eV 2 ). Alternatively it is also possible that v e — > oscillate in the region of 
parameter space suggested by the LSND experiment p2]| . Note however that models with 
3 ordinary and 3 sterile neutrinos may be quantitatively different due to the additional os- 
cillation modes possible. We intend to discuss some of these models in a forthcoming paper 



1 35 



In this four neutrino scenario v T is assumed to be the heaviest neutrino so that v r 
oscillations have 5m 2 s < 0, and create L Vt first at a temperature 



-5m 



7 T« 16 (— v^J MeV - ( 35 ) 

Note that the creation of L Ut implies that the matter term for — > u s oscillations is also 
generated (note that the a^ s {p) term in Eq.(|^) is proportional to LW ~ 2L U + L Ut + L„ e ). 
Thus, the creation of a large L Vr asymmetry also implies the creation of a large function 
which can potentially suppress the — > v s oscillations. For example, for maximal vacuum 
mixing, the matter mixing angle for — > v s oscillations is, 

^ W - l + (a! M - (36) 

In order to simply explain the qualitative behaviour of this oscillation system it is useful 
to consider the quantities (a as ), where (a as ) = a as (p = 3.15T) and p ~ 3.15T is the mean 
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momentum of a neutrino in a Fermi Dirac distribution with zero chemical potential. [Of 
course in our numerical work the momentum distribution will be taken into account]. In 
the region T ~ T c rs where the exponential growth of L Ut occurs, L Ut is generated so that 
(a TS ) ~ 10. If Vp — * v s oscillations do not create significant L v , then 



1 \Sm 2 

IMI * » 1 if \Sm%\ « \5mU- (37) 



Thus, if Up — > v s oscillations do not have time to generate significant L v then the large 
matter term will be generated which will suppress the — > v s oscillations. During the 
subsequent evolution, (a TS ) ~ 1, so that (a MS ) remains very large and thus the — > v s 
oscillations will always be heavily suppressed (until very low temperatures where T <C T/). 
It is important to realise though, that in response to the L Ut created by v T — > v s oscillations 
the U/j, — > u s oscillations can potentially create L v such that — > 0, since this is an 
approximate fixed point of the — > z/ s oscillation system. Note that the evolution of L y 
due to Up — > z/ s oscillations is dominated by the oscillations in the resonance region where 
b^s ~ <Vs (assuming that L Vr > for definiteness)Q. The resonance momentum for maximal 
— > v s oscillations can be obtained from the condition X(p) = 0. From Eqs. (|ll||r9"|) it follows 
that p^ s is given by 

* = (») 

Thus, p^| s is proportional to ~ 2L„ + + L^. As L Vt is created this causes p^ s to 
increase. The — > i/ s oscillations either create L„ fast enough so that p^ s ~ (p) ~ 3.15T 
or they do not. If the — > i/ a oscillations create sufficiently, then (|a MS — 6 MS |) ~ 10. 
Furthermore the ^ — > u s oscillations will continue to keep {(a^ — 6 MS |) ~ 10 for lower 
temperatures because these oscillations become more effective at lower temperatures since 
they are not suppressed so much by the collisions. Thus, we need to compute the evolution 
of L yT , L Ufi due to u T — > v s and — > v s oscillations through the high temperature region 
(-^- ~ T ~ T initia i) in order to obtain the parameter space where v T —>■ v s oscillations 
generate a large L^' which is not subsequently destroyed by — ► 2/ s oscillations. Of course 
the lower temperature (T ~ T c rs /2) evolution of lepton number can also affect BBN in 
several ways, and we will return to this issue in section V. 



As in Ref. [20], we assume that the v T ^v^v a system reduces to two two-flavour oscil- 



lations, v T — > v a and za, — > v s . This simplifying assumption is discussed in some detail in 



Ref. p5| . Heuristically, it is expected that this simplifying assumption is justified because 
the MSW resonance momentum of the two oscillation modes are generally different (and 
different to the v T — > resonance momentum due to the matter effects). Also it should 
be noted that the v T — > u^, v T — > v e and —>■ v e oscillations can be approximately ne- 
glected at high temperatures [Here 'high' means T ~ T C TS which is typically in the range 



10 Of course if L Vt > 0, then the — ► v s MSW resonance occurs for the neutrino oscillations 
(rather than the anti-neutrino oscillations), which create L v < 0. 
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20 ~ T™ /MeV ~ 60], because these oscillations do not create or destroy lepton number 
efficiently at high temperature because there are almost equal numbers of v e ,v^ v T neutri- 
nos. Also v e — > v a oscillations would be expected to be heavily suppressed by the large tau 
lepton number (assuming that |£m 2 J <C \5m% s \)- 

For the v T — > v s and — > v s oscillations we define a density matrix (say P(p) for v T — > v a 
oscillations and Q(p) for — > v s oscillations). Thus, we have 

NuM = \p Q {p) (1 + P-M) N°{p), N Vfl {p) = ±Q (p) (1 + Q z {p)) N°(p), 

N„M = \PoiP) (1 - P-M) N°(p) = \Qo(p) (1 - Q z (p)) N°{p). (39) 



The evolution of the functions Pi(p),Qi(p) are given by Eqs.^^D, where 5m 2 , sin 2 29 = 
5m 2 . s) sin 2 29™ [<5m 2 s , sin 2 29q S ] in Eq.(|TJ) for Pi(p) [Qi(p)] (of course, for the evolution equa- 
tions for Pi(p) [Qiip)], a = t [p] in Eq. (|T9D ). At each time step, the lepton numbers L Ut , L v 
are computed using Eq.(^) [in the case of dL u /dt, Pi(p) — > Qiip)]- The two sets of evo- 
lution equations for Pi(p),Qi(p) are coupled together because they each depend on both 
L Ut and L v . They are also coupled together through the population of the sterile state. 
In particular the population of v s by v T — > v s oscillations can affect —>■ v s oscillations. 
Similarly the population of v s by z/ M — > v s oscillations can affect v T — > v s oscillations. This 
effect can only be important if the population of sterile neutrinos becomes relatively large, 
i.e. n Us ~ 0.1. Note that the population of v s due to z/ M — > v s oscillations is negligible in 
the temperature region T ~ T™. This is because T™ is typically high enough so that the 
—> v s oscillations are suppressed by the matter effects and also damped by collisions so 
that they cannot populate a significant number of sterile states. However, if sin 2 29™ is large 
enough, then the v T — > v s oscillations generate significant and approximately equal numbers 
of v s , u s in the region T ~ T™ . This population of sterile neutrinos will obviously affect the 
Vn — > v s oscillations. The population of v s , v s by v T — > u s oscillations can be incorporated 
by noting that 

N (v) 1 1 

-j^M = 2 P oW(l " PM) = 2^o(p)(l - QMl (40) 

and similarly for the anti- neutrinos. This obviously implies that the rate of change of Qiip) 
due to v T — > v s oscillations (we denote this quantity by §j.Qi{p)\v r -^v s ) satisfies 

| [Qo(p)(i - QM)] * | [Po(p)(i - PM)\ ■ (4i) 

Furthermore, the population of sterile neutrinos by v T — > v s oscillations obviously does not 
directly affect the number of neutrinos. Thus, 

-N u ^p)\ UT ^ Us = - [Q (p)(l + Q z (p))] \»^v. = 0. (42) 
Solving Eqs. ([4~l] ,|4~^) (and the analogous equations for the anti-neutrinos) implies that 
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|o»Mk^ = |q.(p)I_, = i (^(i - p.W) - ^W)) , 

d - d 

= ~2 g (p) " M) — dT oip) ) ' (43) 

where here we have neglected the tiny difference in u s , v s populations (which is actually 
necessary for self consistency). When solving the evolution equations for Pi{p),Qi{p) we 
have included the contribution, Eq. (|43D . 

Our initial conditions for the numerical integration are (as explained in section III), 
p x,y{p) _= Px, y {p) = and P z (p) = P z (p) = P (p) = P (p) = 1 (and similarly for the 
Qi(p), Qi(p), i — x,y, z, 0). We also set the initial values of all of the neutrino asymmetries 
to zero, and took T initia i = 4T™ S . Of course the results do not depend on T initia i so long as 
it is high enough (and T initia i = 4T r T / s is certainly high enough). As we discussed in section 
III (and in Ref. ||20||), the results are also independent of the initial values of the neutrino 
asymmetries so long as they are less than about 10~ 5 . We also utilise the approximation 
of integrating around the region of the MSW resonance [taking / = 7 in Eq.(]2"§D] for both 
u T — > v s and — > v s oscillations. We checked the stability of our results by taking larger 
slices of momentum space. 

Performing the necessary numerical work, we obtained the region of parameter space 
where the created by v T — ► u s oscillations is not destroyed by maximal — * v s oscil- 
lations. This region is given in Figure 2. Of course, we must also require that the sterile 
neutrinos do not become significantly populated by v T — > v s oscillations. Recall from Eq. ([§]) 



that 8N?£ N ~ 0.6 implies the constraint |20[ 



sin 2 29 n ~ 4 x 10" 



eV 2 



\5m 



2 I 

TS I 



(44) 



Note that v T — > v s oscillations populate the sterile neutrinos in the temperature region 
T ~ where a TS (p) is very small. The constraint, Eq.(f44|) is given by the dashed-dotted 
line in Figure 2. There are also other contributions to N^ N , which are due to the low 
temperature evolution of L Vt) L u ^ and L Ue , including the population of sterile neutrinos at 
low temperatures. We will return to these issues in section V. 

Let us now compare our new results with the earlier calculations of Ref. p0| . In Figure 7 
of Ref. |20| we took \5m 2 atmos \ = 10~ 2 eV 2 and computed the allowed region by numerically 



integrating an approximate solution of the quantum kinetic equations. This approximation 
(which we called the 'static approximation' in Ref. p0| ) holds provided that the system 
is smooth on the interaction time scale 1/D(p). For this Sm 2 tmos our new results are in 
agreement with the previous calculation of Ref. pUfl in the region sin 2 2# ~ 3 x 10~ 7 . This 



result was anticipated in Ref. ||20|| , where it was found that the evolution of lepton number 
was smooth enough for the static approximation to be acceptable in this parameter range. In 
the case where sin 26 ~ 3 x 10 our new results show that the allowed region is somewhat 



smaller when compared with the approximation of Ref. |2C . 
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Finally recall that in the context of the standard big bang model, the age of the Universe 
implies an upper bound on the mass of the tau neutrino which is [36] 



< 



100 eV, 



(45) 



m 



and hence 



\8m 2 TS \ £ 10 4 eV 2 



(46) 



Some authors (including myself) have taken a more stringent bound of m VT ~ 40 eV (which 
implies |5m 2 s | ~ 10 3 eV 2 ). While the latter bound certainly appears to be favoured, we feel 
it is prudent to be cautious, and we therefore adopt the weaker bound Eq. (|46|) . 

V Detailed implications for big bang nucleosynthesis 

Let us briefly summarise the story so far. In the region above the solid line(s) in Figure 
2, u T — > v s oscillations create a large L Vt which suppresses — > v s oscillations so that these 
oscillations cannot significantly populate the sterile state (for T ~ 0.5 MeV). Let us denote 
the contribution to NfM from the oscillations, v ^ — > u s [y T — > v s in the region T ~ T C TS ) by 
the notation 5iN^ N (8 2 N^ N ). With this notation, in the parameter region to the right 
(left) of the dashed-dotted line 5 2 N^ N > 0.6 (5 2 N^ N < 0.6). Also, in the region above 
(below) the solid line(s) in Figure 2, S X N^ N < 0.1 (SiN^ N + S 2 N^ N ~ 1). 

It is important to realise that depends on the number densities of u e , v e through 

BBN nuclear reaction rates (as well as the expansion rate). Thus, it is important to study 
the evolution of the four neutrino system down to low temperatures T ~ 0.5 MeV. The 
evolution of this four neutrino system down to low temperatures has already been studied 
in some detail in Ref. plj and we include a discussion here for completeness. Since the 
final value of L Ut , L[ t , is quite big there is a significant modification to the momentum 
distribution of the tau and anti-tau neutrinos. If the tau neutrinos also oscillate into muon 
and electron neutrinos, then some of the tau lepton number can be transferred to the electron 
and muon neutrinos. In other words small L Ue ,L v will be created. The details are quite 
independent of the intergenerational mixing angles so long as they are small (here, small 
means sin 2 29q ~ 0.1 )Q and depend only on 5m 2 . In Ref. |2T| the contribution to N^M N due 



to the modification of the u e , v e momentum distributions (which we denote by d^N^f) was 
found to be 



11 They cannot be arbitrarily small lest the oscillations become non adiabatic (typically they must 
be greater than about 3 x 1CP 10 ). 





0(0), for \5m 2 \/eV 2 ~ 3, 

-0.45(0.45) for 3 ~ \Sm 2 \/eV 2 ~ 1000, 

-0.50(0.50) for \5m 2 \/eV 2 ~ 1000, 



(47) 
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for L Va > (L Ue < 0). Note that 5 3 N B ^ BN in the above equation is actually continuous 
so there is really a smooth transition region between the three Sm 2 regions which we have 
neglected for simplicity^. The creation of the lepton number can also change the expansion 
rate since the energy density of the Universe will be modified a bit by the modification of the 
number and momentum distributions that occur when the large L Vt neutrino asymmetry is 
created. The contribution to NfM N due to this effect (which we denote by 8 4 N B f B ) was 
calculated to be [F2TJ, 



S 4 N BBN ~ for \5m 2 \/eV 2 ~ 3, 

5 4 N B f BN ~ -0.05 for 3 ~ \5m 2 \/eV 2 ~ 1000, 

8 4 N B f B f N ~ 0.40 for \5m 2 \/eV 2 ~ 1000. (48) 

Finally note that the results in Eq.(fI7|) and Eq.(pE8"D are only valid provided that 8\NfM + 

8 2 N B ^ BN ~ 0.1. This is only approximately true in the region above the solid line(s) in 
Figure 2 and in the region where v T — ► v s oscillations do not significantly populate v s at 
high temperature T ~ T c . The region of parameter space where v T — > v s oscillations do not 
significantly populate v s at high temperature (i.e. 5 2 N BBN ~ 0.1) is given by0, 



sin 2 2# n - 1.3 X 10" 



eV 2 
| bra 2 



(49) 



Comparing this equation with Eq.(H) (which assumes that 8iN BB ^ ~ 0.6) we see that it 
is only slightly more stringent. If there is a significant number density of v s coming from 
v T — > u s oscillations at high temperature, then this will reduce the size of L' u . The effects 
given in Eq.fl4"T|) and Eq.(^5|) will therefore be reduced. 

Observe that in the region in Figure 2 below the solid line(s), the sterile neutrino 
will be populated by — > v s oscillations at a temperature T ~ 6 — 10 MeV. When 
this occurs the growth in L Vt is cutoff because the number of tau and sterile neutrinos 
become approximately equal. A these temperatures, \L Vr \ ~ 10~ 2 so that the possi- 
ble modifications of u e , u e are negligible in this case. In other words N B f BN ~ 4, with 

v e// -ru 2 ±y e f f — ±, Weff ' u 4 ly eff 

i , 1 -\ i ■ \- 1 1 1 /". -f if ril m ill in w ■ \- 1 \ 



5 l N BBN + 8 2 N BBN ~ 1, S 3 N BBN , 6 A N$ N « 0.1. 

To summarise, it is possible to identify 4 distinct contributions to N B M N . They are: 



(1) 5tN B f BN . This is the contribution to N B f BN which arises from the change in the expansion 
rate of the Universe due to the population of the sterile neutrinos by z/„ — * u s oscillations. 



12 Actually it should be noted that 5 2 N^ f BN , S 3 N B f ^ N and 5 A N^ f BN are all continuous above the 
solid line(s) in Figure 2. This means that their values change smoothly. The contribution 5\N B j BN , 
on the other hand is discontinuous across the solid line(s). 



13 Note that there is a mistake in the corresponding equation in Ref. [21|. Eq.(|49|) is the correct 
equation. 
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(2) 5 2 Nf f B f N . This is the contribution to Nf f BN which arises from the change in the expansion 
rate due to the population of sterile neutrinos by v T — > v s oscillations for temperatures before 
the exponential growth of the L Vt occurs (i.e. T ~ T C TS ). 

(3) 5-^N BB j N . This is the contribution to N^ BN from the direct modification of the nuclear 
reaction rates (such as v e + iV <-> P + e~, v e + P <-> N + e + ). This occurs because of the 
modification of the momentum distributions of v e , v e due to a small L Ve asymmetry which 
is transferred from L Ut by v T — > v e oscillations. 

(4) S^N B j^ N . This is the contribution to N B j BN from the change in the expansion rate due to 
the modification of the energy densities of the neutrinos at low temperature T ~ T/ s which 
is caused by the large L Ut . 

Note that both the effects (2) and (4) are essentially due to the population of v s by v T — > v s 
oscillations. We label them distinctly because the effect (2) occurs at high temperatures 
before the generation of neutrino asymmetry occurs (in this region u s , v s are populated 
approximately in equal numbers) while the effect (4) occurs because of the population of V s 
(taking L Ut > for definiteness) which necessarily occurs in the temperature region where 
the neutrino asymmetry nears its final value. Finally, the effective neutrino number for BBN 
is related to the contributions, 5iN B j BN , in the obvious way 

C = 3 + D.C (50) 

i=l 

The results of this section are summarised in Figures 3,4. Figure 3 (Figure 4) is for 
the case where L Ve > (L Ue < 0)Q. As these figures show, there is a significant chunk of 
parameter space where NfM < 3. Thus, not only is the — > v s oscillation solution to 
the atmospheric neutrino anomaly consistent with a stringent BBN bound of N B j BN < 3.6, 
but there is a range of parameters where this solution implies that there are less than 3 
neutrinos for BBN. It is interesting that there are some recent indications that N^ BN < 3 
is actually favoured by the data especially if the low deuterium abundance measurements 
D/H ~ 3 x 10~ 5 |J7| in quasar absorption systems are correct P[ ^ ne l esson here 



is that there are well motivated extensions of the standard model which give NfM < 3. 
Thus, if it turns out that N^ N < 3 is required for consistency of the standard big bang 



1 The sign of L Ue is the same as the sign of L Vt , which cannot be predicted at the moment (see 
Ref. [20 1 for some discussion on this point). 



15 It has been argued that N^ BN can be as large as 4.5 [^(| and still be consistent with the 
data. However this paper assumed that D/H could be as large as 2.5 x 10~ 4 which seems to be 



disfavoured by recent studies pj ( see a l so Ref- HI I f° r a review). Nevertheless, things are not 



completely clear at the moment, so one should keep in mind the possibility that N B ^ N as large as 
4.5, as suggested by Ref. [M, may be consistent with the standard big bang model. 
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model, then this should be taken seriously as an indication that sterile neutrinos exist[^]. 

Finally, the scenario we have given in this section holds provided that v e — > v s oscillations 
can be neglected. In the case where |5mgJ <C 1 eV 2 , as happens for the particular case 
studied in Ref. || (where v e — > v s oscillations solve the solar neutrino problem), there is 
no impact for BBN from v e — > v s oscillations. However, in the case where \Sm 2 es \/eV 2 ~ 10 
(which is compatible with the LSND experiment), it turns out that the v e — > v s oscillations 
cannot be neglected for BBN. In particular, consider the possibility that 

m Ve < « m Us < m Vr and 5m 2 es ~ 1 eV 2 , (51) 

along with the requirement that the atmospheric neutrino anomaly is solved by — > v s os- 
cillations, Eq . (|33"D . We will assume that sin 2 29q s is small enough so that v e — > v s oscillations 



are unable to create significant L Ve at high temperatures, T ~ T™. [The situation at high 
temperature is qualitatively similar to the case studied in section IV] . The low temperature 
evolution of L Ue will be dominated by the MSW transitions. It is instructive to consider 
the resonance momentum of v e — > v s oscillations, p e r s es . At low temperatures, where the b(p) 
terms can be neglected, the resonance condition arises from the equation ±a(p) = cos2#o 
(~ 1 in this case for both v T — > v s and v e — > v s oscillations). Taking L Vt > for definiteness, 
which means that at low temperatures, the v T — > v s have a MSW resonance, and the v e — > v s 
oscillations have a MSW resonance (rather than the v e — > v s oscillations). This opposite 
behaviour occurs because 5m 2 s > while 5m 2 s < 0. Thus it is clear that the v e v s 
oscillations generate L Ve with the opposite sign to L Vt . Using Eq.(|l^), 

V f es n 2 5m 2 es 5m 2 es MeV 4 0.23 



0-5-T7f^^T^- (52) 



T 4C(3)v^G F T 4 L( e ) ' eV 2 T 4 m ' 

Eventually, for T ~ T/ r [see Eq.®], L Vr -> V Vt ~ 0.23 (for 3 ~ \5m 2 TS \/eV 2 ~ 1000). For 
these temperatures, w 2L Ue + L[ . Thus, as T — > 0, Eq . (|52D implies that pl s es /T — > oo 
and there will be significant MSW conversion of v e —>■ v s . In this case there cannot be 
complete MSW conversion because of the structure of Eq. (0) . In fact Eq. (|52"D , implies that 
the creation of L Ue actually increases the rate at which p e r s es jT moves through the momentum 
distribution. Thus, L Ve is created much more rapidly than L Vt . Indeed it must be created 
so rapidly that the oscillations are not completely adiabatic which is why incomplete MSW 
conversion occurs in this case. In fact, it is possible to show that as p e r s es /T — > oo, — > 
(approximately), i.e. L Ue — > —Lf /2 ~ —0.12. Since significant generation of L Ue cannot 
occur until temperatures where pl s es /T ~ 1, it follows from Eq . (|52|) that for 5m 2 s <C leV 2 
the significant creation of L Ue occurs at temperatures much less than 0.5 MeV, so that there 



Of course neutrino oscillations is not the only mechanism which can give N^ f B f N < 3, however 
it is perhaps the most well motivated given the existing neutrino anomalies. For some other ways 

rBBl 



to get N B P f N < 3, see for example, Ref. IT 



17 Note that from Eq. (f33[) it follows that and v s are maximal mixtures of two nearly degenerate 
mass eigenstates in this case. 
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is no significant impact for BBN in this case. However, for 5m 2 s ~ 1 eV 2 , the L Ve would 
be created during the time when the N + u e <-> P + e~, P + v e «-> iV + e + reactions are 
still occurring. Thus, in this case there will be a significant modification to BBN. This 
connection between the oscillations in the LSND parameter range and BBN is extremely 
interesting, but we will not present any quantitative numerical results here. However, we 
do intend to discuss this effect more fully in the context of models with 3 ordinary and 3 
sterile neutrinos in the forthcoming paper 



VI Implications for the hot+cold dark matter model 

In the scenario considered in this paper, where za, — > v s oscillations solve the atmospheric 



neutrino anomaly, a BBN bound of < 3.6 implies m Vr ~ 4 eV for |<fr?t^ mos | 

10~ 2,5 eV 2 (see Figure 2). Neutrino masses in the eV range have long been considered to be 
interesting for cosmology. This is because these neutrinos can make a significant contribution 
to the energy density of the Universe. In the standard big bang model, the contribution of 
massive neutrinos to the energy density is given by the well known formula, 

^'"neutrino , n„„ Tr i \""/ 

/r92 eV 

where h is the usual cosmological parameter parameterising the uncertainty in the Hubble 
constant. Thus, neutrinos in the eV mass range are a well known and well motivated 
candidate for hot dark matter. Some studies have suggested that a hot+cold dark matter 



mixture with £l n eutrinos — 0.20 — 0.25 can nicely explain the structure formation ||42|| . If only 
one eV neutrino is assumed, which we take to be the tau neutrino, then the neutrino mass 
favoured in the usual hot+cold dark matter scenario is |42 | 



3 eV ~ m„ T ~ 7 eV. (54) 

This situation is modified somewhat when light sterile neutrinos exist. The reason is that 
the v T — > v s oscillations generate such a large L Vt that the total number of tau neutrinos is 
actually significantly reduced. In the region where 3 ~ \5m 2 TS \/eV 2 ~ 1000, the L{ ~ 0.23 
[from Eq. ([?[)]. The large final lepton number occurs because about 70% of the anti-neutrinos 
have been depleted (assuming L Vt > 0) which means that the total number of tau neutrinos 
(i.e. tau + anti -tau neutrinos) is approximately 0.65 of a standard neutrino species. (Note 
that the total number of neutrinos hasn't changed much, the missing heavy anti-tau neu- 
trinos have just been converted into light sterile states.) This significant depletion of tau 
neutrinos implies that the tau neutrino mass favoured in the hot+cold dark matter scenario 
(which suggests that fl ne utrinos — 0.20 — 0.25) is actually about 50% larger than the naive 
expectation. Thus when v T oscillates with a light u s , the hot+cold dark matter scenario 
actually favours a tau neutrino mass of 

5 eV ~ 77v ~ 10 eV, (55) 

rather than 3 — 7 eV. This means that |<5m^J is expected to be in the range 

25 eV 2 ~ \5m 2 TS \ ~ 100 eV 2 , (56) 
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assuming here that the sterile neutrino is much lighter than the tau neutrino. The hot+cold 
dark matter region, Eq.floq) is the shaded band on Figure 2. From this figure, we see 
that for \Sm 2 tmos \ = 10~ 2 eV 2 there is only a relatively small range of sin 2 29™ where the 
favoured hot+cold dark matter region is compatible with N^ N ~ 3.6. However, for lower 
values of \5m 2 atmos \ which are actually favoured by the superKamiokande data, there is a 
significant range for sin 2 26q s where the favoured hot+cold dark matter region is compatible 
with Nf f B f N ~ 3.6. 



Finally note that detailed studies [|2] of structure formation in the hot+cold dark matter 
model show that they are sensitive to the number of eV neutrinos (not just fl n eutrmos) ■ These 
studies typically assume that the number of eV neutrinos (usually taken to be degenerate) 
at the epoch of matter radiation equality (we will denote this quantity by N^) is an integer, 

jdm 



i.e., Nff) = 1, 2, 3. It is important to understand that this is only true provided that sterile 



neutrinos do not exist. Indeed, as we have just explained above, we expect N^Jj ~ 0.65 in 
the case where only v T is in the eV mass range and oscillates with a light sterile neutrin 



VII Implications for the comic microwave background 

During the next decade or so, high precision measurements of the anisotropy of the 
cosmic microwave background (CMB) will be performed by several experiments (such as 
the PLANK and MAP missions). From these experiments it may be possible to estimate 
accurately the radiation content of the Universe at the epoch of photon - matter decoupling 



H| . In this context it is important to note that sterile neutrinos can leave their 'imprint' 



on the cosmic microwave background |j45| . Thus it is useful to introduce the quantities 
iV^j, Ngff where N^tf (^//) i s the effective number of light neutrinos (heavy neutrinos) 
at the epoch of photon decoupling. In this context, 'light' means much less than about 
an eV (which will make them relativistic at the epoch of photon decoupling) and 'heavy' 
means more than about an eV (which will make them approximately non-relativistic) . Of 
course in the minimal standard model of particle physics which has 3 massless neutrinos, 
^ef B f N = eff = ^' ^eff = 0- However in the four neutrino model case that we are discussing 
in this paper, N? f j N ^ N l e u ff and N^ f ^ in general. [Note that N?fi N ^ N l e u ff not just 
because of the heavy tau neutrino but also because of the contribution of L Ue to N^ N ] . 

If we assume that m„ T ~ 4 eV (as suggested by Figure 2) and m Ue , m v ^m Vs <C leV then 

N % f ~ ^+ N lu ff = P^ + P^ + P-^ (57) 
where (p,) are the mass (energy) densities with n Q (p ) being the mass (energy) density of 



, the 4 neutrino model of Ref. where — > v T oscillations solved the atmospheric 
neutrino anomaly was studied. In that particular model fn, u T are assumed to be approximately 
degenerate and in the eV mass range. In that model the oscillations of the v„ — > v s and v T — * v s 
created a large L M ~ L T ~ 0.16 and Nfj\ ~ 1.5. 



18 In Ref. [24 
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a massless neutrino with Fermi-Dirac distribution with zero chemical potential. Assuming 



that 3 ~ \Sm 2 TS \/eV 2 ~ 1000, usine; results from Ref. [Ell, 

N% f ~ Nfp f ~ 0.65, iV^ ~ 2.2. (58) 

The precise measurements of the CMB may well prove to be quite useful in distinguishing 
between various competing explanations of the neutrino anomalies, since each model should 
leave quite a distinctive imprint on the CMB. 

VIII Concluding remarks 

In this paper we have numerically integrated the quantum kinetic equations to obtain the 
region of parameter space where the — > v s oscillation solution to the atmospheric neutrino 
anomaly is consistent with a BBN bound of 3.6 effective neutrinos. The consistency occurs 
because the v T — *> v s oscillations create a large L Ut asymmetry. This large asymmetry 
modifies the effective potential for — > v s oscillations which become heavily suppressed 
for the range of parameters given in Figure 2. These dynamically induced matter effects 
prevent the z/ M — > v s oscillations from significantly populating the sterile state for the entire 



'allowed region' of Figure 2. This work confirms and improves the previous study [20] which 
utilised an approximate solution of the quantum kinetic equations (which we called the 
'static' approximation in Ref. ||20|1 ). We have also discussed the detailed implications of 
this scenario for BBN (see figures 3,4), the hot+cold dark matter model, and also for the 
forthcoming precision measurements of the cosmic microwave background. 

In this paper we assumed the existence of only one light sterile neutrino. Probably the 
simplest four neutrino scheme which can solve the atmospheric neutrino anomaly by — ► v s 
oscillations and also solve the solar neutrino problem is the case considered in Ref. [[J (where 
small angle MSW enhanced v e — > v s oscillations solve the solar neutrino problem). Our 
results indicate that the tau neutrino mass should be larger than 3 — 4 eV if consistency 
with a stringent BBN bound of 3.6 neutrinos is required. Note however that the scenario of 
Ref. || cannot explain the LSND anomaly f22|. The LSND experiment |22| has provided 



strong evidence for z/ M — > v e and z/ M — > v e oscillations. The suggested parameter range is 
sin 2 28 ~ 10~ 2 and 



0.2 £ \5mi nd \/eV 2 £ 6. (59) 

One may well wonder whether or not there exists any four neutrino scheme which can solve 
all three experimental anomalies (with the — > u s oscillations solving the atmospheric 
neutrino anomaly), and still be consistent with a stringent BBN bound of 3.6 neutrinos? 
The only potential scheme which comes to mind is [PHI H 



m„ ~ m. 



>m v ,m Vs , (60) 



19 Please excuse the sloppy notation in Eq.(60). Of course the mass eigenstates are linear combi- 
nations of weak eigenstates. 
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with v e — > iv oscillations solving the solar neutrino problem, either through vacuum oscil- 
lations or through MSW enhanced oscillations (this means that \5m 2 Te \ ~ 1(T 3 eV 2 ). The 
LSND anomaly is explained by — > z/ e oscillations with 

5m 2 s ~ <5m 2 s ~ ~ 5mf snd ~ 6 el/ 2 . (61) 

In this scenario, v T — > v s oscillations will generate L Vt and u e — > v s oscillations generate 
L Ve . However because the \Smf snd \ is so low, the results in Figure 2 indicate that, except 
for a very small region of parameter space (|5m 2 | ~ eV 2 , sin 2 26 ~ 1O _5 0), we expect 
the ffj, — > v s oscillations to generate a L u asymmetry such that — > 0. This means that 
we expect that the — > v s oscillations will eventually fully populate the sterile neutrino 
at a temperature around 6 — 10 MeV. We conclude that this scheme does not seem to 
be compatible with a BBN bound of 3.6 neutrinos. Hence, if experimental data indicate 
that U/j, — > v s oscillations are required to explain the atmospheric neutrino data, and if the 
solar and LSND anomalies have been correctly interpreted in terms of neutrino oscillations, 
then N^ N < 4 actually suggests the need for more than four neutrinos. Of course, if 
light sterile neutrinos exist then the theoretically most attractive possibility is that there 
are three of them, just like there are three types of ordinary neutrinos and three types of 
charged leptons etc. Furthermore there are well motivated extensions of the standard model 
which have three light sterile neutrinos. For example, parity conserving theories [0] (models 
with similar neutrino phenomenology but without a mirror sector have also been proposed 
in Ref. M) necessarily have three light sterile neutrinos which can naturally explain all 
of the neutrino anomalies if neutrinos have mass (see the second paper of Ref. [0 for a 
summary). These schemes also have the nice property that they can explain these neutrino 
anomalies without assuming any large mixing between generations, and they can also have a 
neutrino mass heirarchy. These last two features are really expected given the situation with 
the quarks and the charged leptons. The implications of these models for early Universe 
cosmology is in progress . 
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Figure Captions 

Figure 1. The creation of \L Vr \/h by v T — ► v s oscillations as a function of temperature, 
T/MeV. In these examples, sin 2 26 = 10" 8 and 5m 2 /eV 2 = -0.5, -50, -5000 for the 
dashed line, solid line and dotted line respectively. 

Figure 2. The 'allowed region' of parameter space in the sin 2 29™, —5m 2 s /eV 2 plane, where 
the v T — > v s oscillations create tau lepton number in such a way as to prevent maximal 
Vn ~> v s oscillations from populating the sterile neutrino for T ~ 0.5 MeV. The bold 
solid line, upper solid line, lower solid line are the boundaries assuming \5m 2 atmos \ / 'eV 2 = 
10~ 2 ' 5 , 10~ 2 , 10~ 3 respectively. The dashed-dotted line is the bound Eq.(|44]) which arises by 
requiring that the v T — > v s oscillations do not populate the sterile neutrino in the temperature 
region before the exponential growth of L Vt occurs. The shaded region is the favoured region 
in hot + cold dark matter models (see the discussion in Section VI). 

Figure 3. Detailed predictions for in the sin 2 2#Q S , — 5m 2 s /eV 2 plane. For simplicity, 

\5m 2 atmos \ = 10~ 2 ' 5 eV 2 is used. In this figure L Ue > is assumed. Note that N^M N is actually 
continuous across the horizontal boundary at 5m? = —1000 eV 2 . The transition region is 
roughly 400 ~ — 5m 2 /eV 2 ~ 3000 (the transition region is not shown for simplicity). 

Figure 4. Same as Figure 3, except that L Ue < is assumed. 
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